The three different phases in the dynamics of chemical reaction networks and their 

relationship to cancer 
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We investigate the catalytic reactions model used in cell modeling. The reaction kinetic is defined 
through the energies of different species of molecules following random independent distribution. The 
related statistical physics model has three phases and these three phases emerged in the dynamics: 
fast dynamics phase, slow dynamic phase and ultra-slow dynamic phase. The phenomenon we found 
is a rather general, does not depend on the details of the model. We assume as a hypothesis that the 
transition between these phases (glassiness degrees) is related to cancer. The imbalance in the rate 
of processes between key aspects of the cell (gene regulation, protein-protein interaction, metabolical 
networks) creates a change in the fine tuning between these key aspects, affects the logics of the cell 
and initiates cancer. It is probable that cancer is a change of phase resulting from increased and 
deregulated metabolic reactions. 

PACS numbers: 87.18.-h, 75.10.Nr 
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I. INTRODUCTION 

When observing the life on the ground of statistical 
physics of complex systems we see a hierarchial level of 
organization and modularity. How can we describe the 
"good, natural" relations between different parts of com- 
plex system? In statistical physics we describe different 
systems through different order parameter, and in equi- 
librium different parts of the system have the same tem- 
perature. For the stability of the system, it is reasonable 
to assume the similarity of order parameter at different 
hierarchy levels and parts (modules) of the system. There 
is a nice similarity with the no-arbitrage condition in fi- 
nancial markets [l| where different stocks can fluctuate in 
equilibrium have only identical order parameters, defines 
as a ratio of driven and diffusion motions. The Random 
Energy Model (REM)[i| and related complexity (replica 
symmetry breaking) order parameters 0] work starting 
from the proteins [J] till quantum chromo- dynamics and 
strings, therefore we assume that the cell organization 
and cancer should not be exclusions, there is a possi- 
bility that REM related ideas can work there. Using the 
known results of REM, we will prove that the phase struc- 
ture of the chemical reaction network kinetics is related 
with the probabilistic distributions of different chemicals 
in the steady state distribution. Then we will speculate 
about the origin of the cancer using similar complexity 
ideas. 

One of the main ideas for cell modeling is to consider 
the dynamical models with network structure [H-Q • One 
can realize this program studying the network of catalytic 
reactions, identifying the phenotype of the cell with the 
attractor of the nonlinear system of differential equations. 



While constructing the kinetic constants, it is advanta- 
geous to have a detailed balance condition. The density 
of different chemical components change according to 
the kinetic constants, and the latter are defined through 
the energies of that components. We are following [9|. 
There are M chemical components. The transformation 
between chemicals Xi and Xj is catalyzed by some com- 
ponent Xc, so has a rate ki^XcXi 

Thus we can write a set of equations 0] 



dxi 



'^W{i,j, c)xciikj^iXi - hjXi) 



(1) 



where W(i,j;c) = W(j,i;c) — 1 when there is a re- 
action and zero otherwise. The kinetic coefficients are 
defined through the energies of the chemical compo- 
nent Ei and some inverse temperature (3 Q kij = 
min{l,exp[P{Ei — Ej)]. Thus the kinetic coefficients are 
defined through the energy landscape. The choice of ki,j 
is quite reasonable for catalytic reaction network: there 
is some finite rate when the reaction goes to the low en- 
ergy configuration, while there is a small probability for 
the reaction in inverse direction. 



The Eq.(l) has a steady state solution 



E,exp[-/?i5,] 



(2) 



In it was considered a homogenous distribution of 
energies in some interval [0, e] and at some temperatures 
glassy behavior was found. The authors claimed that 
the same is the situation in case of normal distribution, 
while the phenomenon lacks in case of log-normal 
distribution or the distribution with the tail. 
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II. RESULTS 



The statistical physics phases of the model. 

Consider a general distribution for energy levels 



1 

2^ 



dh exp[-/i£; + ln(M)0(/i)] (3) 




where 4){h) is a some function. For the normal distribu- 
tion we have = /3^/2. 

When considering pure deterministic dynamics with M 
degrees, we can map our model to the statistical physics 
model REM with In M degrees assuming large M and 
independence of energy distribution. 

In our model Ei are the energies connected with dif- 
ferent chemical components of the model. 

Nevertheless let us introduce a "partition function" Z, 



Z 



M 



The statistical physics model with the 



partition function Z has two phases. At small /? it is 
in the paramagnetic phase (PM) with fast relaxation. 
When the following equation has a real positive solution 



(4) 



the second, spin-glass phase (SG) is possible in the model 
at /3 > /3c- Two phases have different probability distri- 
butions for the Z in the statics . While formulated in 
spin-spin interaction version, in SG phase there is a slow 
relaxation, see also [13 ■ Thus we related some statistical 
physics with the quenched disorder of the dynamic model 

(1)- 

How can we distinguish two phases? At different 
phases there are different distributions functions for Z, 
and also different expressions for the order parameter 

m — 1 ^i(^i)^. 

One has two expressions for the order parameter (con- 
nected with replica symmetry breaking) [isj 



m 



P 



(5) 



For the normal distribution case /3c = Q . While there 
are M degrees in our dynamic model instead of InM 
spins in related REM, we observed the phase transition 
in the dynamics, found first in Q. A similar model has 
been solved in 

In the PM phase the order parameter m{t) converges 
to the steady state value as [I4I 



\m{t) — 1| ^ exp[— fli * t] 
while in the SG phase [l^l 

\m{t) - m| exp[-a2 * t™] 
where ai,a2 are some constants. 



(6) 



(7) 




FIG. 1: The C{t) = ^'^^^l^^^'^" vs time. N = 100. 
Every chemical is connected with K = IQ other chemicals. (a) 
The low line has no ferromagnetic configuration, jS — 0.8/3c, 
the next line has j3 = Pc * 1.4, no ferromagnetic configu- 
rations, the higher lines correspond to the /? = 1.4/3c and 
ferromagnetic configurations with the relative concentration 
0.08, 0.15, 0.30, 0.60. (b)/3 = /3c* 0.8, there is no ferromagnetic 
configurations, the higher lines correspond to the P — 1.4/Jc 
and ferromagnetic configurations with the relative concentra- 
tion 0.08,0.15,0.30,0.60. 
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Let us now add a special energy level Eq = —JqIuM. 
When Jo > /3c, third, ferromagnetic phase (FM), is pos- 
sible 0. We investigated the dynamics of the model in 
this case. When the steady state concentration of the 
"ferromagnetic" chemical is comparable with the total 
concentration of the system, the reaction dynamics is be- 
coming slower than even in the SG phase for the initial 
homogenous distribution. The reason of ultra-slow relax- 
ation is the energy gap between the ferromagnetic level 
and other levels. We performed a numeric for the case of 
normal distribution (realistic for the biology of the cell 
[Hi, [Hi) with M = 100, 7^ = 10 see Fig. 1, and observed 
the phase transition phenomenon in the dynamics. The 
SG phase Fig la is much more sensitive to the appearance 
of FM configuration, than the PM phase. Fig lb. 

This model can be applied for the reactions in the cell. 
We assume that the normal chemical kinetics of the cell 
corresponds to the SG phase, contrary to PM or FM 
phases. 

How general are these results? 

To obtain two phases we assume that: a. The non- 
linear system of equation has one steady state solution 
and b. randomness of steady state concentrations xf. 
The second condition is well established experimentally 
for the concentration of chemicals in the cell [H- The 
concentration of chemicals fluctuate from cell to cell. 
We used a concrete form of the nonlinear dynamic. Let 
us consider another version of nonlinear dynamics, i.e. 
^ = Ej,c^(«'.?''C)(a:c)^(fcj,i2:i - h^^Xi). We performed 
numerics and found the same qualitative picture: the sta- 
tistical physics phase for the quenched disorder defines 
the phase of the dynamic. If we consider the relaxation 
from the initial configuration near the steady state, then 
the results of different nonlinear attractor models cer- 
tainly are similar; they are defined mainly by k{i,j). 

How can we apply these results to the cell-cell inter- 
action, gene expression problem? Ao formulated well the 
idea of adaptive landscape for bio- networks In the 
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Ao's approach, while there is a noise, transversal (mag- 
netic like) field and friction, the system has a given steady 
state defined by a potential landscape. This potential 
landscape should be identified with the fitness. The lat- 
ter sometimes is directly connected with the ordinary free 
ener gy. This idea have been applied to the cancer anal- 
ysis [131 . The robustness and plasticity aspects in the 
dynamics of the cancer cells have been considered in [18 1 
and [l^. The key point here is again the concept of 
the landscape again. Again fast and slow phases in the 
dynamics are possible, for our derivation was important 
just the steady state distribution of energy (fitness). We 
considered SG-PM phase transition at different temper- 
atures. The same transition happens when the steady 
state distributions are changed. This is an important 
observation as the temperature is rather a constant in 
biological systems. 

Connection between the phase in cell's chemical 
kinetics and the type of fitness landscape. 

According to our classification, the steady state could 
be identified either with SG, FM, or PM ones. The other 
characteristic of the model are the basin of attractions 
of these steady states. It is an advantage to have large 
basins of attractions, but simultaneously a slow relax- 
ation dynamics. 

A very interesting point is the connection of behavior 
of the cell (basin of attraction -|- phase of quenched dis- 
order) with the evolution behavior. Here it was found 
that the dynamics of the evolving population drastically 
depends on the character of the fitness landscape func- 
tion: the FM, or SG or their border case, FM-SG 
It is highly intriguing if the character of fitness func- 
tion of evolving cell will be the same as the character of 
the chemical reaction kinetics. A similar ideas about the 
tight relation between cell processes and evolution char- 
acteristics has been suggested in a series of articles by K. 
Kaneko and co-workers, [l^, where they claim that the 
large basin of attraction, related with the given pheno- 
type, can be connected with the mutational robustness. 

In Q also speculations are made about resonance in 
complex systems, assuming identical complexity param- 
eters. It has been assumed that the possibility of para- 
metric amplification of the motion (resonance) and para- 
metric fast attenuation (anti-resonance) is one of the key 
features of normal living systems. Such a property can 
exist for some parameters of dynamic system, without 
being present for other values. 

Speculations about cancer. 

How can we understand the cancer? For a recent re- 
view see (2l| . Here we distinguish two aspects: the origin 
and the dynamic behavior. The evolution dynamics as- 
pect of cancer is well know since [22| - 
of metabolism is also well recognized 
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The key role 
Looking for 

simple and general origin for the cancer, we assume the 
following hypothesis. Since the start of the life, there is 
some fine tuning between different key aspects of the life 
such as gene regulation, metabolism, pH, cytoskeleton. 
It appears that these different aspects of life are inter- 



twined, and their interaction defines the decision making, 
the logics of the cell [1^ . 

The fine tuning of the cell should be understood as an 
existence of proper statistical physics phases, identical 
complexity parameters, oscillations, and this fine tuning 
is supposed to be increasing in parallel with the evolution. 
According to our hypothesis, in the case of cancer: ex- 
ternal aggression like chronic infiammation or increased 
cellular metabolism caused by oncogene activation, dis- 
rupts the degree of the fine tuning, change the reaction 
rates and as a result, the decision making, the logics of 
the cell. 

Having less degree of fine tuning, the cell, in some 
sense, returned back in evolution history, using old ma- 
chinery. A candidate of such a poor fine tuning, can be 
the change in the reaction rates (glassiness), discussed 
in the previous section. The slow, glassy like dynamics 
is typical for the cell to support the non-equilibrium in 
the cell, an important aspect of the life [27|. Another 
side of phenomenon is connected with the memory pro- 
cesses, the glassy (slow) phase possess such memory, 
while the fast phase docs not. The necessity of glassy 
dynamics is well recognized in case of immune system 
[28| . We assume that the cancer cell has another ver- 
sion of glassiness than the healthy cells. It is intriguing 
that the osmotic pressure plays a crucial role for both 
cancer and cytoskeleton glassiness [13]. In case of 
metabolic network the cancer cells choose a fast dynam- 
ics, and perhaps the reaction network system is in the 
PM phase. 

The fast reaction rate and less fine tuning are tightly 
related. Normal cells can either be in anabolic or 
catabolic state. They use oxygen and burn, for example, 
glucose. It results energy, water and Carbonic gas. The 
alternative pathway is anabolism. They use energy water 
and Carbonic gas to synthesize for example glucose. The 
cells can synthesize other compounds like DNA , RNA, 
cholesterol.... In normal cells the two pathways cannot be 
done at the same time. But at different time of the day 
one cell can either be in anabolic or in catabolic phase. 
In cancer cells, some enzymes (modules) are on the an- 
abolic mode, some others are on the catabolic mode. It 
results a deregulation of the metabolism which could be 
seen as a different phase [l^ . 

Consider the change of decision making mechanisms. 
Probably this is the start point of cancer. The metabolic 
rate and decision making are key features of living mat- 
ter [sOIjUH. They should be considered both on a single 
cell leve l, and in a cell-cell ( tissular) interaction [32 1- 
[33| , I34II . We should carefully analyze the metabolic net- 
work [35| and decision making in the healthy cell, and 
compare them with those in the cancer cell. Accord- 
ing to [2^, the cancer cell is less dependent on the sur- 
rounding constraints than the healthy cell. An intriguing 
possibility is the de-synchronization in decision making 
between healthy cells and cancer cells, for example due 
to lack of polarity in cancer cells [s^. Fortunately the 
concept of dc-synchronization is already generalized for 
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non- harmonic processes |37| . 

We mentioned already the hypothesis of Q about a 
property of hving system to reveal both resonance and 
anti-resonance features. In case of cancer the second 
property (apostasis) is lost. In it was suggested to 
construct a rather simplified models of cells so as to 
capture key features of the cancer phenomenon, instead 
of looking into too complicated models. In case we try to 
give a simplified phenomenological model of cancer via 
some dynamical model, we should take care about the 
mentioned feature of the model (a presence of resonance 
and anti- resonance) . 



III. DISCUSSION 

We checked the existence of different phases (degree 
of glassiness) in the chemical reaction network dynamics: 
fast, slow and ultra-slow, connected with the quenched 
disorder. While the first two phases has been found nu- 
merically in , we determined the conditions when they 
emerge and also found the third phase. The phenomenon 
exists only for the special distributions of chemicals with 
a nontrivial solution for Eq.(6), which related to the dis- 
tribution of proteins in the cells [l^. Looking for a 
simple and complexity related key reason for the can- 
cer, we suggested a hypothesis that the healthy and can- 
cer cells have different phases in the chemical reaction 
network dynamics, connected with the glassy properties. 
The glassy like properties of statistical physics models 
arc not artifacts of the modeling, they have very deep 
and universal meaning as complexity order parameters. 
One can apply the spin-glass order parameters, statistical 
physics phases, to the immune system proteins [1], 
cytoskeleton [23| , the reaction network (the current work 
result). In the cells, both normal and malignant, there 
are many modules playing the role of logical units. The 
cancer cell has a poor fine tuning of different modules. 



some of these modules change their meaning, affecting 
the decision making aspects of the cell. We should per- 
form a careful analysis of different aspects in metabolic 
reactions [ssj . gene regulatory network, protein-protein 
interaction network, cytoskeleton dynamics to under- 
stand different choices in case of the healthy cell and 
cancer, the difference in a decision making schemes. We 
should verify our hypothesis, that different aspects of the 
cell (organism) should share the same complexity order 
parameters. We emphasize the idea of "energy" land- 
scape in modern modeling of cell dynamics, including 
population dynamics. We found that the statistics of dis- 
tribution of different degrees at peaks of the landscape is 
as important as the landscape itself. The landscape mod- 
els like the ones should be completed also by a 
description of statistics of different components at peaks. 
We suggest to try experimentally measure both the dis- 
tribution of different chemicals in healthy and cancer cell, 
and also compare the rate of (metabolic) reactions, try- 
ing in this way derive quantitative criteria to identify 
the transition between SG (healthy, slow) and PM (fast) 
phases. Besides identifying different statistical physics 
phases of chemical reaction networks, we suggest to fo- 
cus on decision making synchronization by cells and the 
existence of anti-resonance property. 

It is possible that the investigation of the advanced 
evolution models, incorporating some decision making 
aspects, will improve our understanding of both the can- 
cer and the origin of life. Only gcnetical aspects or 
metabolism are not sufficient to describe the origin of 
life Hi. 
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the support, K. Kaneko, H. Q. Chen, T. J. Kobayashi for 
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